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Abstract 

The formalism that describes radiative-capture reactions at low energies within 
an extended two-cluster potential model is presented. Construction of the operator 
of single-photon emission is based on a generalisation of the Siegert theorem with 
which the amplitude of the electromagnetic process is constructed in an explicitly 
gauge-independent way. While the starting point for this construction is a micro- 
scopic (single-nucleon) current model, the resulting operator of low-energy photon 
emission by a two-cluster system is expressed in terms of macroscopic quantities 
for the clusters and does not depend directly on their intrinsic coordinates and mo- 
menta. The multichannel algebraic scattering (MCAS) approach has been used to 
construct the initial- and final-state wave functions. We present a general expres- 
sion for the scattering wave function obtained from the MCAS T matrix taking into 
account inelastic channels and Coulomb distortion. The developed formalism has 
been tested on the ^He(a, 7)^Be reaction cross section at astrophysical energies. The 
energy dependence of the evaluated cross section and S factor agrees well with that 
extracted from measurement though the calculated quantities slightly overestimate 
data. 

PACS: 21.45.-v; 24.10.Eq; 25.40.Lw; 26.65.+t 

Keywords: Nuclear astrophysics; Radiative capture; Photon-emission operator; Cou- 
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1 Introduction 

In spite of several decades of intensive experimental and theoretical studies of the radiative 
capture of nuclei at low-energies, reactions of this type remain a focus of research activities. 
One of the main reasons for this continued and growing interest in this kind of nuclear 
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processes is their importance in understanding the Big Bang nucleo-synthesis and stellar 
evolution (particularly processes that occur in the interior of the Sun). Another aspect of 
astrophysical relevance that requires improved knowledge of the stellar fusion reactions, 
is the solar neutrino flux problem. In particular, a recent assessment pQ of the current 
understanding of solar neutrino fluxes called for new measurements of the ^He(Q;, 7)''Be 
reaction cross section at very low (solar) energies with an accuracy of better than ±5%. 
Of course, it is a serious challenge for an experiment measuring cross sections at low 
energies to reach such a precision since the reaction yield rate is extremely small due to 
the strong Coulomb repulsion. 

Some current challenges for nuclear astrophysics can be met with new experimen- 
tal facilities such as the ones at the Laboratory for Underground Nuclear Astrophysics 
(LUNA) [.2j that is part of Italy's Gran Sasso National Laboratory (LNGS). The LUNA 
collaboration has reported on their first measurements [3j of the ^He(a,7)''Be total cross 
section at very low (down to 127 keV) center-of-mass (cm.) energies with a total uncer- 
tainty of 4%. High-precision measurements of cross sections for other radiative-capture 
reactions of astrophysical relevance are under way [1]. However, some processes are, and 
will remain the province of purely theoretical studies. Nonetheless, the impressive progress 
achieved in recent times in the experimental studies of low-energy radiative-capture reac- 
tions is an additional argument for boosting research activities on the theoretical side as 
well. 

The photon-emission process relevant to stellar fusion reactions is known [5] to occur 
dominantly at a large (extra-nuclear) distance. In general, then, one may assume that no 
intermediate "compound" nuclear system is formed. If that is so, the process is defined as 
a direct-capture reaction. The ^He(a,7)''Be reaction, which provides termination of the 
proton-proton stellar chain through ^He burning, is a one case where successful application 
of the direct-capture model as a description of the photon emission process has been made. 
One of the earliest detailed analysis of this reaction, and of the mirror process ^H(a, 7)^Li, 
assuming a two-body direct-capture picture was made in Ref. [6_J. In the calculations of 
the total cross section and of the astrophysical 5* factor, the authors of that paper used 
a bound-state radial wave function taken in the asymptotic (Whittaker) form while a 
repulsive nuclear hard-core model that took into account Coulomb distortion was used 
to describe the initial (scattering) state. Such a choice for the initial-state model was 
favoured by the analysis of the then existing experimental data for elastic a scattering 
from the A = 3 nucleus and the phase shifts deduced therefrom. It was shown many years 
later [7| that this repulsive effect can also result from inner oscillation in the overlap with 
the relative wave function of the two nuclei. Those later calculations were made using 
both a two-body potential model and the orthogonality condition model (OCM) [S] to 
construct the wave function. 

The direct nature of many of stellar reactions suggests that using a cluster model 
may be an appropriate approach for their description. In general, a great amount of 
experimental information on these reactions has been accumulated and extensive cluster- 
model studies of the stellar fusion processes have been made [H]. Methodologically, studies 
have ranged from the potential cluster models [TT] ) which consider the two clusters as 
point-like objects in a two-body potential, to microscopic approaches which include the 
resonating-group method (RGM) [12] as well as the generator coordinate method [13] and 
the OCM [g, both of which are closely related to the RGM. 

In addition to the cluster models, the variational Monte Carlo technique has been 
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elaborated [13] and applied in studies of few-nucleon systems. The dynamical input in 
this approach is given by microscopic (nucleon-nucleon) interactions. Doing so provides a 
universality of description and reduces the number of model parameters. In an application 
of this approach [15], total cross sections for the reactions '^He(a,7)^Be and ^H(a,7)''Li 
from calculations were found. But this technique could not reproduce adequately nuclear 
binding energies, and the states of ^Be were so under-bound as to be nonphysical. On 
the other hand, microscopic grounds to justify use of potential cluster models have been 
established within the so-called microscopic potential model j9]. Possible extensions of 
the two-cluster model space were also studied within that approach (see Ref. [16]). 

Nevertheless, many of essential dynamical features of low-energy radiative capture 
processes with light nuclei are still missing in the existing theoretical models. One of the 
important issues that has not been properly addressed so far is the role in the EM process 
of collective (rotational and vibrational) degrees of freedom of the colliding nuclei. This is 
a serious omission since neglect of low-energy collective excitations and the corresponding 
resonances they produce can alter cross sections at low energies substantially. Thus the 
usual theoretical S'-factor extrapolation from the quantities determined by measurement, 
to the energy region of astrophysical interest can be incorrect. 

Recently a multi-channel algebraic scattering (MCAS) approach [17J has been devel- 
oped in which the interaction of a nucleon with a nucleus is described within the collective 
model [IH]. Extension of this approach to the treatment of the two-cluster systems has 
been made [IS]. The MCAS model exploits the so-called orthogonalizing pseudo-potential 
(OPP) method [20] reformulated to take account of the influence of the Pauli exclusion 
principle in the collective model. This method is closely related to the OCM [8J, in which 
the physical two-body states are obtained by renormalization that ensures their orthogo- 
nality to a certain number of unphysical bound states. The OPP method does not involve 
intrinsic degrees of freedom for the two fragments. Instead, an extra nonlocal contribution 
to the original two-body potential is constructed [17]. The MCAS approach solves the 
low energy coupled equations for the interacting systems, and to use information gleaned 
from such in capture reaction studies, the photon-emission operator has to be developed 
within a two-cluster coupled-channel model that also takes into account cluster sizes and 
the explicit coupling to their intrinsic excitations. 

To construct such an operator, we start from the extension [21] of the Siegert the- 
orem [22j, expressing the electromagnetic (EM) current operator in terms of so-called 
generalised electric and magnetic moments of the nuclear system. Then, the amplitude 
of the EM transition is determined by matrix elements of that operator between nuclear 
states and the electric and magnetic field strengths. This explicitly gauge independent rep- 
resentation for the amplitude guarantees the fulfilment of the Siegert theorem for electric 
transitions in the long- wavelength limit. Decomposition of the EM operator into electric 
and magnetic parts is accomplished in a 3-dimensional form which assists in derivation of 
expressions for the low-energy operator of photon emission (or absorption) by a system 
of two extended clusters in terms of their macroscopic properties (the mean square radii, 
magnetic dipole moments and electric quadrupole moments) starting from the microscopic 
(single-particle) treatment of the total system. 

Herein we present expressions for the two-cluster EM operator and for calculation 
of observables for low-energy nuclear radiative capture instigated by this operator. The 
wave functions for the initial (scattering) and final (bound) states will be built within the 
MCAS scheme. In the next section, the manifestly gauge independent expression for the 
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photon-emission amplitude based on the extension of the Siegert theorem is described. 
Sec. |3] contains derivation of the operator of single-photon emission by a system of two 
extended clusters, while the wave function construction according to the MCAS method 
is detailed in Sec. |H Then, in Sec. El results of application of this general formalism to 
describe the ^He(a,7)^Be reaction are reported. Results for the cross section and the S 
factor are discussed and compared with experimental data, and notably that from the 
LUNA experiment [3J. In the context of applications of the MCAS approach, the reaction 
studied is considered as a simplified test calculation since it does not yet involve inelastic 
channels. Our future plans are to use this method to study more complex reactions where 
inclusion of low-lying inelastic channels will be an important aspect of the dynamics. 
Concluding remarks are given in Sec. [61 

2 Photon-emission amplitude based on extension of 
the Siegert theorem 

The amplitude of transition of a nuclear system from an initial state | Pf, i) to a state 
I Pf'^ f) with emission of a single photon with the energy (momentum) E^^ (fc) and the 
polarization vector Syik) = {eo{k), £{k)) is given by 

%f = (27r)'5(P, ~Pf-k) 6{E, - Ej - E^) T,f, (1) 
(z/ = 0,l,2,3) 

where Jjy(O) is the EM current operator at the space-time point x = (t,x) = 0, Ei and 
Pi {Ef and Pf) denote the total energy and momentum for the initial (final) state, and 
indices i and / refer to system "intrinsic" degrees of freedom. 

For the operator Ju{x) and the Hamiltonian H of the nuclear system, the continuity 
equation 

div Jix)=t[H,Joix)] , (3) 
has to hold. This gives rise to the gauge independence (GI) condition 

k{Pf,f\J{0)\P,;t) ={E,-Ef){Pf-f\p{Q)\Pi-i), (4) 

where p(0) = Jo(0), for the transition matrix element. 

In practical calculations, however, this condition is often violated because of short- 
comings in the description of the nuclear dynamics. For instance, for a system consisting 
of Atot nucleons , typically the EM current is given by the one-body contribution [23], 

Atot 

p(o) = = E ^An) , (5) 
€=1 

J(0) = JW = E 5{r,)] + z E [^e X ^in)] , (6) 
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where = e, fi^ = fip (e^ = 0, fi^ = fin) for protons (neutrons), e is the elementary charge, 
fip {fin) denotes the proton (neutron) magnetic moment in nuclear magnetons, and pg, 
and are the momentum, coordinate and spin operators for a nucleon number ^. This 
current does not obey Eq. ([2]), and one instead has 

[Ptot, J(0)] = [K, p(0)] = [H, p(0)] - [V, p(0)] , (7) 

Atot 

-Ptot = 5Z ' 

C = l 

where i^' is the nuclear system kinetic-energy operator and V denotes the nuclear poten- 
tial. Evidently, the GI condition, Eq. (jl]), is not met for the model specified by Eqs. ^ 
and ([2]) if ^ contains nonlocal (momentum-dependent) contributions. An interaction 
current has to be added to restore the gauge independence. 

Furthermore, it may happen that nuclear wave functions used in calculations are not 
eigenstates of H corresponding to eigenvalues Ei and Ef for the initial and final state. 
That is the case in calculation of transition amplitudes for direct radiative capture when 
asymptotic radial wave functions are used for the initial and/or final states. The same 
problem exists in any situation where, due to the increased complexity of the boundary 
conditions, the scattering process is treated with a Hamiltonian that has been simplified 
with respect to that used for the nuclear bound state. 

In view of the difficulties, which are encountered in description of the nuclear system, 
the representation of the amplitude of an EM process discussed in Refs. [2T|, [211 ESI ES] 
can be useful. This representation expresses Tjj in a manifestly gauge- independent way 
in terms of the electric {E{k)) and magnetic {H{k)) field strengths, viz. 

T,f = E{k) D,f{k) + H{k) M.,f{k) , (8) 

E{k) = -i[2{2T;fE.l"^ {{Ei - Ef)e{k) - kso{k)) , (9) 

H{k) = -i[2{2n)^ E^]-^/^k X e{k) , (10) 

where Dif{k) and Mif{k) are matrix elements of the so-called generalised electric and 
magnetic dipole moments of the system determined by the original current J^. This 
relationship provides an extension of the Siegert theorem [22] and is an alternative to 
Eq. ([2]). According to the representation, derivation of which relies upon Eq. ([3]) and the 
diagonality of the total charge operator, we can redefine the current matrix element as 

{PfJ\p{0)\Pf,t) = -tkD,f{k) , (11) 

{Pf, f I J(0) I Pf, i) = -i{E, - Ef) Dif{k) - zM,f{k) X k . (12) 

This provides the automatic fulfilment of the GI condition, Eq. (jl]). 

Decomposition of the EM current operator into the electric and magnetic parts, which 
gives rise to Eq. ([H]), has been discussed in detail in Ref. [21]. Thus we give only a brief 
outline of the corresponding derivation and clarify the meaning of the quantities Dif{k) 
and Mif{k) by discussing their behaviour in the long- wavelength limit. In particular, we 
establish the link between the electric contribution in Eq. ([8]) to the well-known Siegert 
result [221 • 
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So, we factor out the center-of-mass (cm.) motion of the nuclear states, i.e. we use 
the separation 

\Pf,z) =|P.)M), \P^-k■f) =\P,-k)\f), (13) 

which, of course, is always valid at non-relativistic energies. Then, the translation prop- 
erties of the total momentum eigenstates allow us to write 

{Pi-k\= {Pi I e^'=^'°' , (14) 

where i^tot is the cm. coordinate operator. The product of the exponent in Eq. (IT^ and 
the vector e{k) can be decomposed into the "electric" and "magnetic" parts using the 
identity (cf. Ref. [21]) 

^g^fcH.ot ^ r{[p^^^,,£i?^^^e^Afcfi,o.] ^ ^^j^^^ ^ ^ fcje^^'^^-} dA , (15) 
Jo 

where the components of vectors Ptot and Ptot obey the usual commutation rule 

[Ptotj, Rtot,i] = -iSji , (j, / = 1, 2, 3) . 

Analogously, 

^^gifeHtot =eo(^l + zfcPtot £ e'^*^^-' dA) . (16) 

Taking into account the continuity equation, Eq. ([3]), and as the total charge operator is 
diagonal, from Eqs. ([2]), f|T5|) and f|T6|) we deduce Eq. ([8]) with 

A/(fc) = {2nr'{f I D{k) I z) , (17) 
Mif{k) = (27r)-3(/ I M{k) I t) . (18) 

The operators 

D{k) = ^\P^ - Xk I Ptot [p(0), H] I Pi) dA , (19) 

hi — hf Jo 

Mik) = -(27r)3 f\Pi - Xk I Ptot X J(0) I Pi) XdX (20) 
Jo 

act in the space of the intrinsic variables of the nuclear states. It should be stressed that 
Eq. (IH]) is quite general and can be used to describe various EM processes with nuclei at 
low and intermediate energies. In case of the inverse reaction (absorption of the photon 
with momentum k), one has to replace Xk by —Xk in Eqs. (1191) and (!20|) . The results 
given by Eqs. ([2]) and ([8]) are identical if the GI condition of Eq. (jlj) holds, but can be 
different otherwise. We choose the representation of the amplitude given by Eq. ([8]) since 
it provides automatically the GI of calculations and helps to reproduce predictions of 
the low-energy theorems based on the gauge invariance. Furthermore, as will be shown in 
Sec. [3|, this form of the amplitude is a convenient starting point for construction of effective 
EM interactions for a system of extended nuclear objects (clusters) with an account for 
their internal structure. 
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The representation through Eq. ([8]) was successfully employed in our studies [27] of 
pion photo- and electroproduction off nucleons and light nuclei. The current paper ad- 
dresses construction of EM operators for applications at very low energies (such as nuclear 
radiative-capture reactions at energies of astrophysical interest). 

Operators defined by Eqs. (HM and (|2UI) act in the space of states characterised by 
intrinsic coordinates 

r'^ = r^- Rtot (21) 



and momenta 



Ko, , (22) 



with commutation relations 



^tot 

€=1 



771 

K,,p'J=^(5,.-^)5,, (23) 



(■C) — 1)2,..., y4tot ; J) ^ — 1) 2, 3) . 

To clarify the meaning of the quantities specified in Eqs. f|T9|) and fl20l) . consider the 
current model as defined by Eqs. (jS]) and After substitution of Eqs. and into 
Eqs. f fT^ and fl2U]) and taking into account relations given in Eqs. fl?T]) - (1251) . we obtain 



D{k) = ^ f'dX{[D,^,{Xk),H,, 
hi — hf Jo 



V ^0 

\k(2P,-\k)^ ,,,, Pi-\k ,,,,, 
+ ^ ^Ant(Afc) + ^^J^Plnt{Xk)} , (24) 

Atot 

Pint(Afc) = Ee€e~*"""'S (25) 

D,^,{Xk) = Y^e^r'^e-'^'^^'^ (26) 

for the electric generalised dipole operator and 

M(fc) = M°'^{k) + M^P''^(fc) , (27) 

rl 2P- — AA* 

M-^fc) = XdX{M-^{Xk) - X Ant(Afc)} , (28) 

M«P^'^(fc) = £ e^s^e-''="'« + k / AMA £ e^(r'5S5)e-'^'="'^ , (29) 

(A/c) = E -^{/'e,e-^^^^'4 , (30) 
l'^ = r'^ X p'^ (31) 



for the magnetic one. Obviously, the second term in Eq. (1291) does not contribute to the 
amplitude in Eq. ([8]). Moreover, for a fc-congruent frame [Pi x k = 0) which will be 
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implied below, contributions from the second term of Eq. and from the integrand in 
Eq. (128|) cancel when emission of a real photon is considered. In that case one usually 
uses the Coulomb gauge, 

Eoik) = , ke{k) = , (32) 

for which the last term in Eq. ( 12^ also vanishes. Partial cancellation between the two 
terms of Eq. ([8]) reflects the fact that, in general, decomposition of an EM transition 
amplitude into electric and magnetic parts is not unique (see, e.g., discussion in Ref. [25]). 
Thus, result, we have 

D{k) = ^ / A„t(Afc) dA , (33) 

At, 



Ei-Ef Jo 



M{k) = / M^^iXk) XdX + Y: e^s^e-'^^'^ , (34) 
Jo 

where i?-"* (-E'™*) is the intrinsic energy of the initial (final) nuclear system. Calculation of 
the amplitude of interest is reduced to calculation of the matrix elements of operators -Dint? 
M°^l^ and M^^^^ between the intrinsic states of the nuclear system. From Eq. (!33l) the 
classical Siegert result ^22j for the electric transitions in the long-wavelength limit A; — >■ 
is readily deduced. In this limit, the operator M{k) becomes that of the total magnetic 
moment for the system. Unlike the common procedure [23l [28] of obtaining the electric 
and magnetic multipole operators, Eq. ([S]) with quantities Dif and Mif determined by 
Eqs. fl33p and decomposes the amplitude into the electric and magnetic parts without 
using partial-wave expansion, and no approximations valid only for low photon energies 
are made. The operators given in Eqs. fl33|) and fl34|) are three-dimensional involving an 
integration over A that can be easily carried out. In the next section we show that this 
representation in terms of generalised electric and magnetic dipole operators, is of partic- 
ular use in derivation of photon-emission operators for cluster-type nuclear structures. In 
fact, with Eqs. (133!) and general expressions for the operator of photon emission by 
a system of two (or more) finite-size composite objects can be derived while taking their 
internal cluster structure into account. 



3 Operator of single-photon emission by a two- 
cluster system 

Assume now that the nuclear system consists of two subsystems (clusters) A {A nucleons) 
and B {B=Atot—A nucleons) and is described by a wave function determined by subsystem 
intrinsic states | 0^) and | 0^) and a vector | \1/ab) corresponding to the di-cluster relative 
motion. We can introduce for such a system a new set of coordinates and momenta as 
follows 

P'.^P^ + ^P. P'.-P^^^. (36) 
n^R.-Rs. P^MIZJ!^, (37) 
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(« = l,2,...,/l; /3 = 1,2,...,S) 

where -Ra,b (-Pa,b) and Ma,b are the coordinate (momentum) and mass of the subsystem 
A, B; and r^, (p^, p^) refers to the "intrinsic" coordinate (momentum) of a nucleon, 
which belongs to cluster A, B. 

For the two-cluster relative coordinate R and momentum P and the subsystem in- 
trinsic operators, the following commutation relations hold 

[Rj, Pi] = iSji , [rff,ptf] = - T^) Sji , 

[rtf^ pit] = Hf, Pi] = [R^,ptf] = , (38) 
(^, K=1,2,...,A [B] for cluster A [B]; j, / = 1, 2, 3) . 

They can be obtained straightforwardly from Eqs. fl35|) - fl37|) with the help of Eq. fl23|) 
and manifest the reciprocal independence of the intrinsic coordinates and momenta for 
the two different subsystems from each other. 

Then, substitution of Eqs. ( l35i) and ( 136|) into Eq. ( l26l) gives 

Ant(Afc) = dA(Afc) e-^^^'=^ + dB(Afe) e^^^'=^+r>AB(Afc) , (39) 
dA,B(Afc) = I: e^rf •^e-^^'^^''' , (40) 

i^AB(Afc) = ^pA(Afc) e--^'^« - ^pB(Afc) R e-^'^^ , (41) 

PA,B(Afc)=|:e^e-'^'»^''' . (42) 

It should be stressed that Eq. ([33]) for the generalised electric dipole operator, with 
iI>int(Afc) determined by Eqs. fl5^ - fH2]) . is valid to exactly the same extent as is the 
one-body current of Eqs. ([5]) and (j6]). Eq. (!39l) provides a convenient starting point to 
construct electric transition operators of different form depending on a particular sit- 
uation. Consider, for instance, low-energy photon emission by a system of two charged 
clusters with substantially different sizes. In addition, due to a strong Coulomb repulsion, 
on average the emission (absorption) process occurs at a very large separation distance, 
considerably greater than the characteristic size of the bigger cluster. Then, one can make 
a low-energy expansion of Eq. fl39|) with respect to specific dynamical variables, that is 
independent of the remaining variables. In other words, it is possible to obtain evaluations 
at different perturbative orders with respect to the three separated expansion parameters 
kr^, kvp and kR. 

In general, the quantity <iA,B(Afc) in Eq. (!39l) gives rise to a generalised electric dipole 
operator for the subsystem A, B and contains all multipole contributions responsible for 
the electric transitions in this subsystem. Analogously, the operator £)AB(Afc) provides 
(for all orders in A;) the electric transition contribution due to the relative motion of the 
two clusters. But, as follows from Eq. 0391) . the separation between the "intrinsic" (for 
the subsystems A and B) transitions and the ones due to the relative di-cluster motion is 
not complete. The "intrinsic" electric operators d^iXk) and dB(Afc) in this formula are 



9 



coupled to the exponential factors dependent on the two-cluster relative coordinate, while 
the quantity £>ab(AA;) includes cluster form-factor operators pxi^k) and p^i^Xk). 

However, such a separation is exact in the long-wavelength limit in which, from 
Eqs. (ESD and ([31]), one gets 



D{0) = + + D^^ , (43) 



= cCeiR , (44) 

^ MbZa-MaZb ,,,, 
" M ' ^^^^ 



where d^'^ = dA,B(0) is the usual electric dipole operator of the system A, B 
Za (-^b) is the component A (B) charge, and D^^ is the electric dipole operator for two 
point-like clusters. This proves that the constructed electric operator fulfils the Siegert 
theorem |i22j at — 0. Note also that there is no recoil in this limit. Therefore the energy 
factor in Eq. (1551) is equal to 1. 

If the states | xa) and | xb) possess a definite parity, and the intrinsic parities of 
the clusters do not change in the EM transitions, by neglecting terms of order O(fc^) in 
Eq. fl39|) we obtain 

pint pint ^ ■ ^ 

D{k) = ^-r^ {d''' - lD^\k)) , (46) 

Df{k) = Y. h {t% + tl + eC^,i;f ) , (48) 
I 

tr = f.h{^4:?-tf-hVt'']') , (49) 



T^"" = SRjRi - 5,iR^ , (50) 

MgZA + MjZB 
Ce2 = j^, , (51) 

(j,/ = 1,2,3) . 

Here, we have taken into account the gauge condition, Eq. ( !32l) . The operator in the 
parentheses of Eq. f HHl) is the electric quadrupole tensor operator for the two-cluster 
system. 

Of course, one can proceed with analysis of Eq. (!39|) and obtain higher-order contri- 
butions to the electric operator defined by Eq. (!33|) . In particular, a correction to the 
two-cluster electric dipole operator D^^ emerges due to the finite cluster sizes. This 
correction can be of great importance when the radiative capture occurs with different 
symmetric nuclei, namely, those having equal numbers of protons and neutrons. For 
those cases, the quantity given by Eq. fHSl) is vanishingly small. Examples are the radia- 
tive capture of ^He by ^^C, by ^^O, or by some heavier nuclei. That capture process is 
very important in assessing stellar evolution and nucleo-synthesis. For simplicity, assume 
that the reaction is analysed within a potential cluster model with no coupling to intrinsic 
excitations of the subsystems A and B (i.e., they are considered as "frozen"), and that the 
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cluster electric quadrupole moments are zero. Then, for sufficiently small photon energy 
{kr^'^ -C 1), from Eq. (!39|) the correction 

Cei ^ CEi{k) = Cei + PACei , (52) 

^^^^ = "18 MZTm ' ^ ^ 

is found. Therein {r\ g) is the mean square charge radius of the cluster A, B. However, the 
corresponding contribution to the operator D{k) has to be taken into account together 
with other terms in Eq. fl33|) of the same order in k. 

The same reasoning can be also applied to the magnetic operator given in Eq. 
Substituting Eqs. ([35D and fl5^ into Eq. fl5U]) gives 

M°f (Afc) = m°^^(Afc) e-'^^''^ + m°^^{Xk) e'^^''^ + M^^{Xk) 



{cZA(Afc) X V, e"'-' M ^{dB(Afc) X V, e*^ m 



2M' ^ ' ' 2M 



R X 'i;A(Afc)] e"*^ M [R x VB{Xk)] e'^^"^ , (54) 



M ' ^ " M 



A,B 



m 



A,B 



(Afc) = ^ 7r^{Zf e~'^'="^^} , (55) 



:^ 2m^ 



- 2Mm!^^(^^) {X,e---} + ^^PB(Afc) {X,e---} ,(56) 

^A,B(Afc) = ^Eeat.f^e-^^'^^e'''} , (57) 
^ €=1 

^« = ^ . ^ = inr- ' = ,^ , 58 

^ iwAB JM 

if'^ = r^A,B ^ pA,B ^ L = RxP. (59) 

The structure of the ffist line in Eq. (15^ is very similar to the structure of Eq. (!39|) 
for the electric transition operator. However, the magnetic operator couples the cluster 
"intrinsic" EM transitions to the relative motion between clusters in a more complicated 
way than does the corresponding electric operator. That is so because it also contains 
vector products of the subsystem intrinsic operators <iA,B and i^a.b with the two-cluster 
relative momentum and coordinate. 

Similarly, the magnetic spin operator of Eq. fl29l) cast in terms of the two-cluster 
coordinates defined by Eq. fl35l) becomes 

M^P^"(fe) = mf "(fc) e-^^*^^ + mf''{k) e'^'=^ , (60) 
mf^ik) = Y.e^s,e-^''^t\ (61) 
where the Coulomb gauge, given by Eq. (!32|) . has been assumed. 
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As was done in the analysis of the operator D{k), we assume that the magnetic 
operator defined in Eq. (l34l) acts in a space of intrinsic states possessing a definite parity. 
Also consider the part of this operator which gives rise to transitions that do not change 
parity. Then, in the limit k ^ 0, all but the first three terms in Eq. can be omitted 
and we can write 



M(0) = fi^ + At' 



AB 



.A,B 



C 



1 

5=1 



M 



1 



Ml 



- /^nC'mi-C' 5 



(62) 
(63) 

(64) 
(65) 



Act \Ma Mb 

where /xn = e/2mN is the nuclear magneton and mN is the nucleon mass. Therefore, 
the generalised magnetic dipole moment determined by Eqs. flMl) and fl^U]) in the long- 
wavelength limit is the total magnetic dipole moment operator for the two-body system. 
This is a sum of the subsystem intrinsic magnetic dipole operators fi^'^ defined in the 
usual way [23l [28] with an additional contribution M^^ that is due to cluster-relative 
orbital motion. Evidently, in the specific case of a two-cluster potential picture without 
coupling to "intrinsic" excitations for the subsystems, transition matrix elements of op- 
erators fi^'^ are determined by the static magnetic dipole moments of the clusters A and 
B. 

We can proceed with the analysis of the magnetic operator given in Eq. flM|) for 
low photon energies in a similar manner to that we have used for the electric operator. 
Disregarding terms of order Olfi-^k"^), the part of this operator which does not change the 
cluster intrinsic parity can be written as 



M{k) = M^'^' - iM^^\k) 
M^" = M(0) , 



(66) 
(67) 



M^'''(k 



+ -Rx 

6 

1 

+ 3 



= {kR) 

Mb 



M 



M, 



B A 



Ma 



fJ-""] + fXNCM2{L,{kR)} 



M 

HA,df\k) 



M 
Ma 
M 



'^dfik) - ^df{k)) X V , 

^A^B; = ^ (^kitf{^ + ekjZ A,B rl^B 



C 



M2 



3AotM \Ma 



M| 
Mb 



(68) 

(69) 
(70) 



where [Hb) and r\ (r^) are the cluster A (B) intrinsic Hamiltonian and the operator 
of the square charge radius, respectively. Details are given in Appendix |Al If we consider 
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the case of a cluster-like potential model, with no coupling to the intrinsic excitations and 
with vanishing static electric quadrupole moments, Eq. (!68|) reduces to 



M*«(fc) = (fcR) + f,„CM2{£,(*fi)} 

^l_mz.(rl) -M,Z,(rl) ^ x V . (71) 

Finally, note that the formalism outlined here paves the way for the development of an 
extended cluster model for description of low-energy photo-nuclear processes taking into 
account the internal structure and EM excitation of the clusters. Starting from the mi- 
croscopic (single-particle) model defined by Eqs. ([5]) and ([6]) for the EM current operator, 
we have derived the low-energy EM transition operators, and specifically separated the 
terms acting on the part of the wave function that describes the relative motion between 
clusters. 

Conversely, the cluster intrinsic structure is described by operators of cluster magnetic 
dipole and electric quadrupole moments and squares of charge radii. Usually, these in- 
trinsic operators are defined microscopically, but the evaluation of their matrix elements 
can be done within any model, be it microscopic (e.g., by the shell model, or by few- 
body techniques for light nuclei) or macroscopic (a collective-type model). Therefore, our 
approach can be considered as a step towards construction of a unified framework for 
treatment of static electromagnetic properties of nuclei and low-energy nuclear radiative 
capture with a link (as shown below) to the description of nuclear two-body scattering 
within the multichannel approach of Ref. [17]. Furthermore, this formalism can be easily 
generalised to the case of three-cluster and more complex systems through the proper 
extension of the set of the Jacobi coordinates and momenta. 



4 MCAS wave functions for two-cluster states 

Consider the radiative-capture process 

A + B ^ C + 7 . (72) 

Initially there is a scattering state | ^l/p''; Ja^a; 0jg</B"^B) characterised by the two- 
body cm. momentum P and a nucleus A (B) intrinsic state (^Jb) with the total 
spin Ja (Jb) and its projection mx {jti-q). The process forms a bound state | ^'j^Jc'^c) 
which has total spin Jq and projection niQ. Assume that the nuclear system is described 
within a two-cluster picture without any direct reference to single-particle coordinates 
and momenta. This conforms to the current MCAS approach [17], which utilises the OPP 
method [20l [29| to account for the Pauli principle. With partial-wave decomposition, the 
two-cluster initial- and final-state wave functions can be written as 

{R I ^p^;0j^Ja'^a;0Jj3"^b"^b) = {a",m'i,m'^ \ 3^jmj(^r,^p) | a,mA,mB) 

J,a,a", 
m" m" 

X I <l>j,,J'lml) I <l>\Jl^m'^) ^'J:\R) , (73) 



13 



and 

{R\<!>%Jcmc) = E C{L'M'^J'^m'j,J'M'Ji,m'^Jcmc) 



a',M',M'^, 



X I <I>\J'a^'a) I <I>IJWb) ^Jca'iR) YL'M'.i^R) , (74) 

wherein we have used the short notations 

\a) = \P; ((L, Ja) J'", Jb) J), 
\a') = \P'; J',J^) Jc) , 

\a")^\P"; iiL",J'i)J",J^) J) , (75) 
for the channel states. We have also made definitions 

(a", ml, m'^ \ yjuA^R^ ^p) I ^a, ^b) = ^ C{L" M'lJ'lmlJ" M" J^m'{,JMj) 

M,M" ,M"', 

xC{LMLJAmAJ"'M"'JBmBJMj) Yl.m'i{^r) Ylj^j^i^p) , (76) 

C(jl"^lj2"^2j3"^3j4m4j5"^5) = (jl?Tlij2?^2 | Js^^s) (j3"^3j4"^4 | js^^s) • (77) 

These formulae are given for a general case when non-central nuclear forces [30] and in- 
elastic {P" 7^ P) two-body channels are to be taken into account. The angular-momentum 
coupling scheme implied is shown explicitly in Eqs. fl75p . 

The coordinate-space scattering radial wave function \I^^„^~*'-* (i?) has an asymptotic 
behaviour 




<"i''\R) R^\jlj^Yt^ \l^,Si"a0^uhP"R) - 5a"aO[-\PR) \ , (7J 

where 5*;^,/^ is an element of the partial-wave S-matrix and O^"-* and O*^^ are the incoming 
and outgoing Coulomb waves respectively. For the latter, we imply [32] 
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0[^\PR) ^ ^^[PR~vH2PR)-^f] ^ (79) 

r] = riiP) = ZaZ^Mab/P , (80) 
with a being the fine-structure constant. Next, we can write 

0'l\pR) =Ol^-\pR) e^'^^(^) , 

0^j\){P"R) =0'!}^\P"R) e-'-L"{P") ^ (81) 

0'1^^\x) = Gl{x)±iFl{x) , (82) 
where, as usual, for the Coulomb regular {Fl) and irregular (Gl) functions, one has 

Ltt 

Fl{PR) ^ sm[PR -r^ln{2PR) - — + a^iP) ] , 

Gl{PR) cos[PR - T] \n{2PR) - ^ + aL{P) ] . (83) 
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The Coulomb phase shift depends on the Sommerfeld parameter t] as 



ai(P) = argr(L + l + ir/(P) 



Using the definitions given in Eq. (IHTl) and introducing the "reduced" S matrix Sf [30] 
through the definition 

^/„„ = e--(^")5,V„e-(^), (84) 
the asymptotic condition fl78|) becomes 



2i 



^5) 



The wave function \[^^„j,^^(i?), which is regular at i? = and satisfies the asymptotic 
condition, Eq. fl85l) . can be constructed from the MCAS multi-channel scattering matrix 
given in the Coulomb-state representation [33] for a separable two-body potential of rank 
N, namely 



N 



Va"aiP,q) = J2iP\ Xa"n)K^{Xan \ Q) , (86) 
n=l 

in which the momentum-space potential form factors are the Fourier-Coulomb transforms, 

' " (87) 



Xanip) = {P I Xan) = \ / F^ipr) Xan(r)dr , 

V 7rp Jo 

of the corresponding coordinate-space form factors. 

In our present study, we employ the procedure [M] to construct the scattering radial 
wave function for the separable potential (l86l) in coordinate space. That involves a direct 
resolution of the coupled-channel Lippmann-Schwinger equation in a finite-rank matrix 
form using the Green's function expressed through the Coulomb functions and for the 
physical outgoing solution, i.e., that obeying condition given in Eq. ( l85i) . Following this 
procedure with the MCAS scattering matrix [17], we get 



FL{PR)6a"a - TtMabv/W $il+)(i?) 



N 



n,n'=l 



+ 0^!i'\P"R) Xa"n{P") } [a - gl^'Ynn' ' ^^^^ 



where 



/ 2 1 roD 

X^jP,P) = i/-- / GL{Pr) Xan{r)dr 

\ H P JR 

xLiP,R) = y-p PLiPr) Xan{r)dr ^ 

Xan{P)^XaniP,0) ■ 



(90) 
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The quantity in the square brackets of Eq. (189!) is a matrix, with its elements in the space 
of separable-decomposition indices being 



(+) 



2M, 



AB 



open »oo 

E 

L a=l 
closed 



^-^^^ P^-x^ + .0 ^-'^^^ 



- J2 Xan{x)^::^—y^Xan'{x) dx 

n -^0 



where 



(91) 



(92) 



P = ^2MAB(^c.m. - ^*'^) , /ia = pMj^B{E^^ - E,,^) , (93) 

where E^^ is the threshold energy of channel a and Ec.m. is the elastic-channel kinetic 
energy. Both are considered in the cm. frame. The function given in Eq. fl89|) oscillates 
at large radii with amplitude given by the multi-channel T matrix. Thus we can write 



S 



o 



R 

J, a" a 



6„ 



(94) 
(95) 



where the asymptotic outgoing Coulomb wave Of^f^'*(oo) is given by Eqs. (1791) and (IHTl) . 

Thus, we have constructed the radial wave function for the two-fragment initial (scat- 
tering) state, which is general in the sense that possibilities of coupling to inelastic channels 
are taken into consideration along with coupling between elastic channels due to the non- 
central part of the two-body potential. The wave function construction for bound states 
can be fulfilled through solving the coupled-channel Schrodinger problem analogously to 
the procedure outlined above for the scattering wave function. That is so as the MCAS 
approach utilises sturmian expansions to construct the separable interactions, with which 
determination of both the positive-energy (resonance) and negative-energy (bound) two- 
body states solutions of the Lippmann-Schwinger equations can be found based upon 
analytic properties of the sturmian functions. As all details have been given [TT], we 
do not reproduce them here. We conclude this section by noting that the inversion of 
the matrix A — Qq^^ contains all the dynamical features of the two-cluster interaction 
process, including possible formation of a compound system at the given resonant energy. 



5 Reaction ^He(a,7)''Be 

The formalism outlined above encompasses studies of EM transitions in which the role of 
the finite sizes and collective excitations of the fragments can be taken into account. How- 
ever, for an initial application of the formalism, we have chosen the reaction '^He(a, 7) ''Be. 
In this case there are no low-energy inelastic two-body channels, whereas the elastic chan- 
nels are uncoupled. Such a choice has been governed in part by the amount of experimental 
information available as well as the many theoretical studies made of it in the past. Fur- 
thermore, a goal of this illustrative calculation was to establish the degree of applicability 
of the MCAS model in a description of EM processes. We also wish to compare theoretical 
ingredients of this calculation with those given by other approaches. In the first instance 
we have not sought to achieve a quantitative agreement with experimental data. Thus 
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we have not varied parameters of the two-body potential in our calculation. Instead, the 
radial wave function's for the initial and final states were constructed using the nuclear 
potential of Ref. \X9\ which had been found to best reproduce the states and resonance 
widths of ^Li. 

The differential cross section for the radiative-capture process can be written as 
d^^-"^- - 4vr^a ^^J^ (| A, ^ - | kD^j T + I M., x k ^ 



dn^ \ P\M + E^ 



eD,f = {f\D{k)\i), eM,j = {f\M{k)\t), k 



+2Re[D*f [Mif X fcjj j , (96) 
A; 



k 



All the quantities are in the cm. frame and averaging (summation) over initial (final) 
nuclear system spin states is implied. Details for calculation of the matrix elements 
{f\D{k)\i) and {f\M{k)\i) are given in Appendix [Bl When calculating the magnetic 
transition matrix element, only the first two terms of the M2 operator, Eq. fl68l) . are 
retained. The last term of that formula is negligibly small for the energy range covered 
by the calculation while the third one is zero for the reaction considered. 

To calculate the radial overlap integrals for the ^He(a,7)''Be reaction, scattering- 
and bound-state radial wave functions have been found following the MCAS formalism 
described in Sec. |H To minimize any impact of computational uncertainties in evaluation 
of the bound-state wave functions upon the results, we replaced these wave functions at 
a very long two-cluster distance (| |> 10 fm) by the properly renormalised Whittaker 
functions that correspond to experimentally known binding energies. The radial wave 
function for the ^Be ground state gives rise to the rms charge radius of 2.61 fm, which is 
in fair agreement with the experimental value of 2.52 ± 0.3 fm [35] • The Coulomb wave 
functions were calculated using the cernlib program [36] . 

The main contribution to the ^He(a,7)^Be reaction amplitude at low energies is well 
known to result from the electric dipole transition. This is so because the quantum 
numbers of both the '^Be ground state (| ) and its excited (| , 429 keV) state allow 
the electric dipole transition to occur from the scattering s-wave. On the other hand, 
the operator in Eq. p7|) responsible for this transition in this case is neither vanishing 
nor small. As a consequence, the angular distribution obtained is rather flat at very low 
two-cluster relative energies. Some small non-uniformity is caused by the interference 
between the dipole transitions from initial s and d waves. At higher energies, the E2 
transition due to the operator given in Eq. (1501) also plays a role, making the cm. cross 
section asymmetric about 90°. Magnetic transitions arising from the operator in Eq. fl66|l 
are only minor contributions to the cross section as shown in Fig. [H These observations 
agree with the conclusions given in Ref. [7J. 

The real parts of the radial overlap functions (Eq. flll6p in Appendix [B]) for the dom- 
inant electric dipole transition for the captures to the ^Be ground state from the relative 
s-state at cm. energies of 100 keV and 1 MeV are shown in Fig. [21 For convenience, 
the Coulomb-phase exponential factor in Eq. (!88|) has been omitted. The initial- and 
final-state radial wave functions constructed within the MCAS scheme exhibit oscilla- 
tions in the radial overlap at small (< 4 fm) distances. These wave functions and the 
resulting overlap function are in remarkable agreement with those given by the OCM and 
potential- model calculations [7j. In Fig. [2] we also compare those results with the overlaps 
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0.5 1.0 1.5 
Ecm. (MeV) 

Figure 1: Multipole contributions to the total cross section for ^He-'^He capture leading 
to the ground state of ^Be. 

when the initial-state radial wave function is a pure Coulomb function and the scatter- 
ing wave function is that for the repulsive hard-core (RHC) model [6]. The hard sphere 
radius was taken as 2.8 fm. These overlap functions indicate that the photon emission 
occurs, on average, at very large extra-nuclear distances. However, as seen also in Fig. [2l 
taking into account the nuclear forces leads to a significant repulsive effect which is very 
important for quantitative estimates. Comparison of the curves for -Ec.m. =1 MeV and 
100 keV shows the degree of spread of the EM interaction area with the decrease in the 
two-cluster relative energy. 

In Figs. [3] and H] respectively we show our results for the total cross section (Jtot and 
for the corresponding astrophysical S-factor which is defined as 

S{E,,^,) = Ec.m.c^tot(^c.m.) exp [2nr]{P)] . (97) 

Here, 

P = v'2MAB^c.m. • 

Our calculated results reproduce the energy dependence for the ^He(Q;, 7) ''Be cross section 
and for the S'-factor very well. However, the calculated cross sections are ~40% larger 
than the data. This is typical [TJ [12] of both the potential and microscopic RGM cluster- 
model calculated results. How this overestimate depends upon the choice of the nuclear 
potential has been studied by Kajino [15]. He discovered that the RGM calculation with a 
particular nuclear force, viz., the MHN interaction ^46], could reproduce the absolute data 
for the '^He(a, 7)^Be cross section and astrophysical 5* factor that were known at that time. 
Nevertheless, the S'-factor extrapolation to the zero energy S{0) = 0.50 ± 0.03 keV-barn 
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Figure 2: The real part of the radial overlap function for the electric dipole transition to 
the ^Be ground state from an initial relative s-wave. The thick and thin solid curves were 
obtained with the scattering radial wave function defined in Eq. flSSl) for cm. energies 
of 1 MeV and 100 keV, respectively. The dotted and dot-dashed lines are the overlaps 
found using the repulsive hard-core (pure Coulomb) initial-state wave functions for a cm. 
energy of 1 MeV. 
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Figure 3: Total cross section of ^He(a,7)''Be reaction. The filled circles at -Ec.m. = 127, 
148 and 169 keV are recent LUNA data [3]. Other experimental points were taken from 
Refs. [37l[38l[39lll0lliTlll2lll3llll]. 
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Figure 4: Astrophysical S factor for reaction ^He(a;, 7) ''Be. The notation is as used in 
Fig. [31 The dotted curve displays the result obtained with the initial-state radial wave 
function constructed for the RHC model. 
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suggested in Ref. |15] is smaller than 5'(0) = 0.547±0.017 keV-barn given by the latest R- 
matrix fit to all current data, including that from the LUNA experiment [3]. Evidently, 
the two-cluster a+'^He configuration does not exhaust all the possibilities for the ''Be 
bound state [IH]. In principle, that could explain why the magnitude of the calculated 
cross section is larger than experiment. However, it was indicated [TB] that extension of 
the two-cluster model space for ''Be by inclusion of the p+^Li channel might even make 
the agreement between the theory and experiment worse. 

In contrast, the variational Monte Carlo calculation [15] made with the two-body 
wave functions constructed from nucleon-nucleon interactions found an S factor that lay 
significantly below the experimental values for the ^He(a,7)^Be reaction. In particular, 
the value of S{0) = 0.40 keV-barn has been extrapolated. But this approach, without 
the inclusion of multi-nucleon forces, was not able to reproduce satisfactorily the ''Be 
bound-state levels and the binding energies for the separated ^He and "^He fragments. 
In addition, phenomenological procedures had to be used to construct the two-cluster 
(7-body) scattering wave functions. With respect to this last point, note that recent 
few-body techniques have been developed to deal with this difficult problem [17] . 

Though our treatment is not of the "ab initio" type, both the initial- and final-state 
radial wave functions have been constructed using the same two-body potential [19] ; that 
potential with which the ^Li two-cluster binding energy and a number of resonance posi- 
tions and widths were determined. This potential also gave the a-|-^He system resonance 
parameters listed in Table 1 and they, too, are in good agreement with experiment. In 
this table, the state positions are given relative to a — ^He breakup threshold and all 
units are MeV. Nevertheless, the ^Be | excited state appears to be too weakly bound 



Table 1: Spectrum of states for ^Be as a + ^He system. 



J" 


Level position (width) 


Exp. level position (width) 


2 


-1.53 


-1.59 


1- 
2 


-0.84 


-1.16 


7- 
2 


3.07 (0.180) 


2.98 (0.175) 


5- 
2 


5.09 (1.2) 


5.14 (1.2) 



in our calculations. This affects the corresponding bound-state wave function behaviour 
in the EM interaction region. As a result, we obtained the branching ratio (ratio of cross 
section for capture into the | state to that for transition into the ground state) equal to 
0.5. That is slightly above the average experimental value of ~0.45. Our calculated ratio 
essentially is independent of the energy in complete agreement with both experimental 
data and other calculations. 

Using the MCAS two-cluster potential [19] the existing data on '^H(a, a)^H and 
^He(^He, ^He)^He scattering reactions for the cm. energy range between 2 and 7 MeV is 
well reproduced. It also gives good representation of the low lying states in the mass-7 
compound nuclei that can be formed. When choosing the "optimal" potential parameters 
in that treatment, there were no evaluation of the nuclear scattering phase shifts at lower 
energies. Comparison of our result for the astrophysical S factor with that obtained in 
case of replacement of the MCAS initial-state radial wave function by the RHC scattering 
wave function (see Fig. H]) suggests that the agreement of our calculations with experimen- 
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tal data could be somewhat improved by correction of the MCAS potential parameters 
for better description of nuclear scattering at very low (below the lowest resonance) ener- 
gies. Though, one should keep in mind, when performing this comparison, that the RHC 
initial-state radial wave function is not properly normalized. 

Adjustment of the potential given in Ref. [T^] to find a better fit to the '^He(a, 7)^Be ra- 
diative capture data has not been attempted in this study as our primary goal has been to 
present a comprehensive formalism for the treatment of radiative-capture reactions within 
the MCAS two-cluster picture. We have made a first apphcation of such an approach, but 
only in a simplified test case. The general approach has been developed, however, with 
a view of application to radiative-capture processes involving more complex nuclei. For 
this reason we paid attention to particular aspects of the EM transition such as the role 
of finite fragment sizes and collective (rotation and vibration) excitations. Nevertheless, 
with the MCAS scheme and without any extra adjustment of the given potential model 
we could reproduce the ^He(a,7)''Be cross section as well as other theoretical models, be 
they microscopic or based on a potential specifically constructed to describe the capture 
cross section. 

6 Conclusions and outlook 

A formalism has been presented to describe direct radiative-capture reactions at low ener- 
gies within an extended two-cluster potential model. Unlike customary potential cluster 
models in which the interacting nuclear fragments are usually treated as point-like objects, 
the EM current operator has been used so that finite sizes of clusters and their intrinsic 
excitations are taken into account. Inherent in the approach is the explicit gauge inde- 
pendence and reproduction of the low-energy behaviour of the EM transition amplitude, 
which follows from the current conservation viz., fulfilment of the Siegert theorem [22] 
for electric transitions in the long- wavelength limit. These properties are ensured by 
extension [2n [25] of the Siegert theorem which served as a starting point to build the 
photon-emission operator. The formalism was used to construct the electric and mag- 
netic operators responsible for low-energy EM transitions in which the intrinsic parities 
of clusters remain unchanged. Starting with a microscopic single-nucleon current model, 
construction of these operators lead to expressions in terms of macroscopic quantities 
which refer to fragment EM static properties or intrinsic transitions such as the cluster 
mean-square radius and magnetic dipole and electric quadrupole moments. Such a fea- 
ture allows one to perform calculations with wave functions of a nuclear model including 
collective-type coupled-channel dynamics. 

Indeed, it has been shown how the photon-emission operator can be utilised to calcu- 
late the cross section of the low-energy radiative-capture process (!72|) using multi-channel 
radial wave functions for the initial (scattering) and final (bound) states. One of the goals 
of this study was to match the constructed EM current with the MCAS coupled-channel 
approach [TTJ US] developed for description of nuclear scattering. We have presented 
the explicit construction for the scattering radial wave function based on the MCAS for- 
malism. The corresponding expression is given for a separable two-body potential of an 
arbitrary finite rank in terms of the potential form factors, the MCAS T matrix and the 
regular and irregular Coulomb functions. This wave function construction is general in the 
sense that it is given for a two-cluster nuclear potential with non-central forces. Possible 
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inelastic channels are taken into account as is Coulomb distortion. Detailed formulae for 
the matrix elements of the EM operator between the constructed nuclear states have also 
been given. Thus, applicability of the MCAS theory has been extended to cover nuclear 
radiative-capture reactions including ones of astrophysical relevance. 

Finally, we have applied the formalism to calculate the total cross section and astro- 
physical S factor for the ^He(a,7)''Be radiative capture at low energies. This reaction is 
strongly dominated by the electric dipole transition, with the effects due to cluster sizes 
and intrinsic structure being negligible. In addition, there are no inelastic channels in 
the a + '^He system, while the elastic channels are uncoupled. Thus many features of the 
developed formalism are not required in analysing this reaction. However, not only does 
this application serve our purpose, the reaction itself is important from the astrophysical 
point of view. For the latter reason it has been extensively studied both theoretically 
(in different models) and experimentally. Our calculation aimed at establishing the level 
at which the MCAS two-body potential scheme is able to reproduce the observed cross- 
section (and ^-factor) energy dependence at low energy without adjustment of potential 
parameters. To obtain the ^He(a,7)''Be cross section, we used the MCAS potential [12] 
with parameters providing the best fit to state positions and resonance widths for the mir- 
ror two-fragment system a -|-^H. The energy dependence of our calculated cross section is 
in a fair agreement with experiment and with that found using other theoretical models. 
The calculated S-factor is about 40 % above the experimental data, a result which is also 
typical of calculations performed within other models, in particular those based upon the 
RGM. 

As shown in Ref. |15], better agreement with experiment can be achieved if one em- 
ploys a two-body potential with parameters tuned to describe fairly well not only the 
two-fragment bound-state levels and resonance parameters, but also the static EM prop- 
erties of the final nucleus. However, the corresponding potential analysis should also 
include the correction to the nucleus static magnetic moment due to the interaction cur- 
rent contribution. That comes from nonlocalities (dependence on the velocity) of the 
two-cluster potential. Besides the "one-body" EM operator discussed in the present pa- 
per, a procedure for constructing the interaction current for the nuclear two-fragment 
system is needed. That is work in progress. 

The formalism presented herein, based on the construction of the EM operator and 
the wave functions built using the MCAS approach [TTl [19] , establishes a single potential 
framework for a unified description of low-energy nuclear scattering, radiative-capture 
reactions, and of static EM properties of nuclei, in the case that they can be treated 
as two-cluster systems. It is our intention to consider scattering and capture reactions 
involving more complex nuclei than the relatively simple '^'^He. In those future studies, 
we expect that the role of cluster sizes and intrinsic excitations in forming the energy 
dependence for the corresponding astrophysical S factors will be of great importance. 

The authors express their gratitude to K. Amos, for critical reading of the manuscript. 
L.G.L. would like to thank the University of Padova and I. N.F.N, (section Padova) for 
the kind hospitality and support in 2006 - 2007. 
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A Appendix: Details of the magnetic operator 

Assume that nuclear system dynamics is described by a Hamiltonian 

H = Ha + Hb + Hab, (98) 

with Ha (Hb) being the subsystem A (B) intrinsic Hamiltonian, and Hab describing the 
two-cluster relative motion. Let us obtain the fc-order correction to the leading-order 
contribution (!62l) in the part of the magnetic operator which gives rise to transitions 
without changing the parity of the cluster intrinsic states. 

First, consider the joint contribution of order fi^k given by the first and sixth terms 
of Eq. flMl) , for which we have 



Xi6(Afc) 



-iX 



a=l 



If the intrinsic potential for the cluster A is local, then the relationship 



HA.rt 



(99) 



(100) 



is equivalent to the one given for the quantities and by Eq. (158|) . However, Eq. fllOOp 
can be considered as a definition of operator which enters Eq. f l54p in a more general 
case (see Ref. [23], p. 393). 

Using Eq. fllOOl) . we get the relationship 



i[HA,dA{Xk)] 



A 
a=l 



Ha, r^e-'^'^"'^ 



-i\kr„ 



HA,e 



-iXkr^ 



a=l a=l 



On the other hand, we can write for the quantities (Ia^b 

cZA,B(Afc) = dA,B(0) - adi%(fc) + 0{e 



(101) 

(102) 
(103) 



(J,^ = l,2,3) 



where t^{^ is the cluster A, B electric quadrupole operator defined by Eq. 



and 



A,B 



' A,B 



A,B 



(104) 



is the square charge radius operator for the corresponding subsystem. Due to the parity 
selection, the contributions from t'A(O) and <iA(0) vanish for transitions between same- 
parity Ha eigenstates. Hence, owing to Eqs. fllOip and fll02p . we can write 



E ^o.vt {kr, 



a=l 



Y^e^rt {kvt)+i[HA,Sl\k) 
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(105) 



omitting higher-order contributions. The latter relationship allows us to write 



a=l 



a=l 



a=l 



"9 E ^ 

^ a=l 



+ 2 



HA,d1\k) 



(106) 



Substituting Eq. (11061) into Eq. (199!) and omitting longitudinal terms i.e. those propor- 
tional to vector k, we get 



Xi6(Afc) = -.^A^(fci2) mr(0) + ^A^ilx [ifA,d«(fc) 



(107) 



Analogously, the joint contribution from the second and seventh terms of Eq. flM|) is given 
by 



X27(Xk)=t-X^(kR) m°'''(0) - -X^Rx 



H^.dg\k) 



Since transitions between states of the same parity cannot be made with the operators 
^a,b(0), the 4**^ and 5**^ terms in Eq. (1541) of the lowest (first) order in g, provide the 
contribution 



V (M.^ ., M^Sl\k)-Mp^4\k) 
A45(Afe) = iX — X V 



(109) 



Substitution of Eqs. ( I107P - 01O9P and the same-order contribution given by Eqs. (!56l) 
into Eq. (IMl) and taking into account the spin contribution coming from Eq. ( l60l) gives 



iM''\k) = -i{kR) 



V M ^ 



Ma b 



+ -i? X 

6 



M 



HA,Sl\k) 



Ma 
M 



H^,4\k) 



. M3dy(.)-MA<(.) ^ ^ _ ^ e 
3 M Im^ 



(110) 



1 



'M, 



M2 



SAtotM V Ma 



M 



'iir 



B Appendix: EM transition matrix elements 



In this Appendix, we present the detailed formulae for the transition matrix elements of 
the EM operator constructed in Sec. [3] between the initial (scattering) and final (bound) 
states described in Sec. IHfor the radiative capture process, Eq. fl72]) . The notation is as 
used in those sections. 

Spherical components of the matrix elements of the operators given in Eqs. (l46l) and 
(166|) between the states specified in Eqs. (173|) and (ITill can be written as 



^int ( 



if 



+ 



if 



E2 



AB 



if 



;ii2) 
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(/ I M^{k) I i) 



4 



if 



+ 



/if 



-i(-l)"A;_ 



Mb 
M 



if 

Ma 



if 



fifR^\^^ + fi^CM2[{h,R.}ij) ' (113) 



if M 

(e, /€ = -!, 0,1) 

where, for simplicity, we ignore the last two terms of expression fl55]) . which either do not 
contribute at all, or are negligibly small for the energy range covered by our calculation 
of this amplitude presented in Sec. [51 Then, the matrix elements 



which enter Eqs. (11121) and (11131) . can be written as 



C{LOJAmAJAmAJB'mBJ,mA + itib) J L [O^J^,^'^ 



(114) 



;il5) 



a,a'a",J 



where 
and 



3 



2j + l 



1 if A is projectile, and B is target 

(—1)"^ otherwise . 



Here, the angular-momentum quantization axis is taken to be in the direction of the 
momentum P. Denoting the radial overlap integral as 



^a'a"a 



for the angular- momentum coupling scheme defined by Eq. flTSl) . we have 



(116) 



L <iia'a"a '^A^A '^b^b^ ' 

x(J,mA + mB,l,^ I Jc'^c) I J/ ^/ I 



J" J 



Jc 



J' 



;ii7) 



X 



(J,mA + mB,2,^ + K I Jcmc) I "^"^ 2 } | 



J" J'^ J 
Jr. 2 J' 



x(0);||QA||0^^)4ti°^ (118) 



JJcmc 
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JJcmc 
a' a" a 



L" 


J'a 










J' 


2 


J) 


{ Jc 


2 





x{L"020\L'0)l'fJ,'-^\ 



(120) 



JJcmc 



X 





L' 


J") 




J'b 






1 




1 Jc 


1 





(121) 



JJcmc 



5rj"h'L"5j'^Jl{-lf '^^''^~^'~'''''^ (J,mA + mB,l,e I ^c^nc) 



X 



J' J 

Jc, 1 



(122) 



X (J, mA + me, 1, e I ^cmc) | f | | f f 'j, ] I'fJ^ , (123) 



II 111 Til IT" 

JcJ{l) 



B J' J" L" (L"010 I L'O) 



JJcmc 
a' a" a 

x(0^, II Ma II 0^^) /„t;:i'^E^ (1^1/^ I 9,i + 

9 

{ J" J' J \ \ ^ ^ ^ 

x{J,mA + mB,g,C + K \ Jcmc) { j ^ , } { L' J'j^ J' 

[Jc 9 J \\l>i ji^^ r 



(124) 



JJcmc 



5j/^j.(-l)-^ +^ +-^A J/ jn L" (L"010 I L'O) 
L" 

^ J' 1 L' 

X ^(-1)^^ (l^lK \ g,^ + K,){J,mA + mB,g,^ + K,\ Jcmc) < 



r 1 1 

Jc 

'II 
B 



J" J 



(125) 



[{^^^Jiair = <^.;.^<5.^,.^^(-i)'^+'^'+''"+^"+''^+^Bjv"L"(L"oio I mi'fj,^:^ 



L" 


J'a 










J' 


9 






9 


'A 



X (|(-1)^L' ^U{U + 1) I |, I +L" v/^"(^" + 1) 



1 1 9 
V L" L" 



,(126) 



where quantities 1 1 M.a,b \ \ 4>^'i^ ) and 1 1 Qa,b | | ) are the reduced matrix 

elements of the magnetic dipole and electric quadrupole operators for cluster A, B. In a 
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particular case of transitions that do not change the fragment intrinsic states (0^/'^ = 

), they are related to the static magnetic dipole (/ia b) and electric quadrupole (Qa b) 
moments of the nucleus A, B through the corresponding definitions 



(>^A,B>^A,B 1 I '^A.B'^A.b) / iA,B II .V II ,A,B 



J' 



Ma,B 1 1 (J) J' ) = AtN/^A,B 



A,B 



^ J' -"'iA,B \ \ Y J' 

■'a.b ■'a.b 



(^A,B^A,B 2 I ^a.B-^A.b) / ,A,B m ^ m ,A,B \ ^ n n'7\ 

■ ■ — {(pj! ^ II Qa,b II <Pj! J = eQA,B ■ (127) 



A,B 



^A,B "'A.B 
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